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Reheating- volume measure for random-walk inflation 
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The recently proposed "reheating-volume" (RV) measure promises to solve the long-standing prob- 
lem of extracting probabilistic predictions from cosmological "multiverse" scenarios involving eternal 
inflation. I give a detailed description of the new measure and its applications to generic models of 
eternal inflation of random-walk type. For those models I derive a general formula for RV-regulated 
probability distributions that is suitable for numerical computations. I show that the results of the 
RV cutoff in random- walk type models are always gauge-invariant and independent of the initial con- 
ditions at the beginning of inflation. In a toy model where equal-time cutoffs lead to the "youngness 
paradox," the RV cutoff yields unbiased results that are distinct from previously proposed measures. 



I. INTRODUCTION AND MOTIVATION 

It was realized in recent years that in many cosmolog- 
ical scenarios the fundamental theory does not predict 
with certainty the values of observable cosmological pa- 
rameters, such as the effective cosmological constant and 
the masses of elementary particles. This is the case for 
the "landscape of string theory" [l|, 0, Q (see also the "re- 
cycling universe" and for models of inflation driven 
by a scalar field (see e.g. for early work). A com- 
mon feature of these cosmological models is the presence 
of eternal inflation, i.e. the absence of a global end to 
inflation in the entire spacetime (see Refs. 0, H, 0] for 
reviews). Eternal inflation gives rise to infinitely many 
causally disconnected regions of the spacetime where the 
cosmological observables may have significantly differ- 
ent values. Hence the program outlined in the early 
works [lol [Til . [T^ was to obtain the probability distribu- 
tion of the cosmological parameters as measured by an 
observer randomly located in the spacetime. The main 
diffuculty in obtaining such probability distributions is 
due to the infinite volume of regions where an observer 
may be located. 

Since the spacetime during inflationary evolution is 
cold and empty, observers may appear only after reheat- 
ing. The standard cosmology after reheating is tightly 
constrained by current experimental knowledge. Hence, 
the average number of observers produced in any freshly- 
reheated spatial domain is a function of cosmological pa- 
rameters in that domain. Calculating that function is, 
in principle, a well-defined astrophysical problem that 
does not involve any infinities. Therefore we focus on 
the problem of obtaining the probability distribution of 
cosmological observables at reheating. 

The set of all spacetime points where reheating takes 
place is a spacelike three-dimensional hypersurface [H, 
[isl . [l^ called the "reheating surface." The hallmark fea- 
ture of eternal inflation is that a finite, initially inflating 
spatial 3-volume typically gives rise to a reheating sur- 
face having an infinite 3-volume (see Fig. [T]). The ge- 
ometry and topology of the reheating surface is quite 
complicated. For instance, the reheating surface con- 
tains infinitely many future-directed spikes around never- 



Figure 1: A l-|-l-dimensional slice of the spacetime in an eter- 
nally inflating universe (numerical simulation in Ref. p2|). 
Shades of different color represent different regions where re- 
heating took place. The reheating surface is the line separat- 
ing the white (inflating) domain and the shaded domains. 



thermalizing comoving worldlines called "eternally inflat- 
ing geodesies" [H, O [l3 ■ It is known that the set of 
"spikes" has a well-defined fractal dimension that can be 
computed in the stochastic approach [l^. Since the re- 
heating surface is a highly inhomogeneous, noncompact 
3-manifold without any symmetries, a "random location" 
on such a surface is mathematically ill-defined. This fea- 
ture of eternal inflation is at the root of the technical and 
conceptual difficulties known collectively as the "measure 
problem" (see Refs. [Sl, |9j, ll8|, ll9|, l20i |23| for reviews). 

To visualize the measure problem, it is convenient to 
consider an initial inflating spacelike region S of hori- 
zon size (an "//-region") and the portion R = R{S) of 
the reheating surface that corresponds to the comoving 
future of S. If the 3-volume of R were finite, the volume- 
weighted average of any observable quantity Q at reheat- 
ing would be defined simply by averaging Q over R, 



(1) 
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where 7 is the induced metric on the 3-surface R. 
This would have been the natural prescription for the 
observer-based average of Q; all higher moments of the 
distribution of Q, such as (Q^), {Q^)j etc., would have 
been well-defined as well. However, in the presence of 
eternal inflation^ the 3-volumc of R is infinite with a 
nonzero probability X((/>o), where = 0o is the initial 
value of the inflaton field at S. The function has 
been computed in slow-roll inflationary models [l5| where 
typically X{<j)Q) « 1 for 0o not too close to reheating. In 
other words, the volume of i? is infinite with a probability 
close to 1. In that case, the straightforward average (H]) 
of a fluctuating quantity Q{x) over R is mathematically 
undefined since y^dPx = 00 and /jj Qy/^cPx — 00. 

The average (Q) can be computed only after imposing 
a volume cutoff on the reheating surface, making its vol- 
ume finite in a controlled way. What has become known 
in cosmology as the "measure problem" is the difhculty 
of coming up with a physically motivated cutoff prescrip- 
tion (informally called a "measure") that makes volume 
averages (Q) well-defined. 

Volume cutoffs are usually implemented by restrict- 
ing the infinite reheating domain i? to a large but finite 
subdomain having a volume V. Then one defines the 
"regularized" distribution p{Q\V) of an observable Q by 
gathering statistics about the values of Q over the finite 
volume V. More precisely, p{Q\V)VdQ is the 3- volume of 
regions (within the finite domain V) where the observ- 
able Q has values in the interval [Q,Q + dQ]. The final 
probability distribution p{Q) is then defined as 



piQ) = lim piQ\V), 

V— s-oo 



(2) 



provided that the limit exists. 

Several cutoffs have been proposed in the literature, 
differing in the choice of the compact subset V and in 
the way V approaches infinity. It has been found early 
on (e.g. 0, that probability distributions, such as 
p{Q), depend sensitively on the choice of the cutoff. This 
is the root of the measure problem. Since a "natural" 
mathematically consistent definition of the measure is ab- 
sent, one judges a cutoff prescription viable if its predic- 
tions are not obviously pathological. Possible pathologies 
include the dependence on choice of spacetime coordi- 
nates [2^ [25! , the "youngness paradox" jglj^ and the 
"Boltzmann brain" problem [11 [H, M, MMM, H HI . 

The presently viable cutoff proposals fall into two 
rough classes that may be designated as "worldline- 
based" and "volume-based" measures (a more fine-grained 
classification of measure proposals can be found in 
Refs. [13, [3) ■ The "worldline" or the "holographic" 



measure [36|, |37[ avoids considering the infinite total 3- 
volume of the reheating surface in the entire spacetime. 
Instead it focuses only on the reheated 3-volume of one 
i/-region surrounding a single randomly chosen comov- 
ing worldline. This measure, by construction, is sensi- 
tive to the initial conditions at the location where the 
worldline starts and is essentially equivalent to perform- 
ing calculations with the comoving-volume probability 
distribution. Proponents of the "holographic" measure 
have argued that the infinite reheating surface cannot be 
considered because the spacetime beyond one iJ-region 
is not adequately described by semiclassical gravity [37| . 
However, the semiclassical approximation was recently 
shown to be valid in a large class of inflationary mod- 
els [38] . In my view, an attempt to count the total volume 
of the reheating surface corresponds more closely to the 
goal of obtaining the probability distribution of observ- 
ables in the entire universe, as measured by a "typical" 
observer (see Refs. [34, 39, 40, 41, 42] for recent dis- 
cussions of "typicality" and accompanying issues). The 
sensitive dependence of "holographic" proposals on the 
conditions at the beginning of inflation also appears to 
be undesirable. Volume-based proposals are insensitive 
to the initial conditions because the 3-volume of the uni- 
verse is, in a certain well-defined sense, dominated by re- 
gions that spent a long time in the inflationary regime.^ 
Existing volume-based proposals include the equal-time 
cutoff d, 1^, [13] , the "spherical cutoff" [l^l , the "comov- 
ing horizon cutoff" [ij, |4^, |4^, the "stationary mea- 
sure" (21. . i47l . the "no-boundary" measure with volume 
weightin g |48| . l49l . [sol . [Slj . the "pseudo-comoving" mea- 
sure [3ll. l52l|. and the most recently proposed "reheating- 
volume" (RV) measure [sst . 

The focus of this article is a more detailed study of 
the RV measure in the context of random-walk eternal 
inflation. As a typical generic model I choose a scenario 
where inflation is driven by the potential 1^(0) of a min- 
imally coupled scalar field </). In this model, there ex- 
ists a range of <j) where large quantum fluctuations dom- 
inate over the deterministic slow-roll evolution, which 
gives rise to eternal self-reproduction of inflationary do- 
mains. I extensively use the stochastic approach to infla- 
tion, which is based on the Fokker-Planck or "diffusion" 
equations (see Ref. 9] for a pedagogical review). The 
results can be straightforwardly generalized to multiple- 
field or non-slow-roll models are straightforward since the 
Fokker-Planck formalism is already developed in those 
contexts [s^ [sj. Applications of the RV measure to 
"landscape" scenarios will be considered elsewhere. 

An attractive feature of the RV measure is that its con- 
struction lacks extraneous geometric elements that could 



Various equivalent conditions for the presence of eternal inflation 
were examined in more detail in Refs. [l5l I23II and [l4j . Here I 
adopt the condition that X{(f)) is nonzero for all <j) in the inflating 
range. 



^ It has been noted that 3-volume is a coordinate-dependent quan- 
tity, and hence statements involving 3-volume need to be formu- 
lated with care |43|| . Indeed there exist time foliations where the 
3-volume of inflationary space does not grow with time. The 
issues of coordinate dependence were analyzed in Ref. |23|| . 
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introduce a bias. An example of a biased measure is the 
equal-time cutoff where one considers the subdomain of 
the reheating surface to the past of a hypersurface of 
fixed proper time t = tc, subsequently letting tc ^ oo. 
It is well known that the volume-weighted distribution of 
observables within a hypersurface of equal proper time is 
strongly dominated by regions where inflation ended very 
recently. A time delay 6t in the onset of reheating due to a 
rare quantum fluctuation is overwhelmingly rewarded by 
an additional volume expansion factor cx exp [SiJinax*^^] j 
where iJmax is roughly the highest Hubble rate accessi- 
ble to the inflaton. This is the essence of the so-called 
youngness paradox that seems unavoidable in an equal- 
time cutoff (see Refs. [Ill and (s^l for recent discussions). 

Moreover, the results of the equal-time cutoff are sen- 
sitive to the choice of the time coordinate ("time gauge") . 
For instance, the proper time can be replaced by the fam- 
ily of time gauges labeled by a constant a, [2J| 



t 



(a) 



(3) 



which interpolate between the proper time (a ~ 0, i(o) = 
t) and the e-folding time (a = 1, = Ina). It has been 
shown that the results of the equal-time cutoff depend 
sensitively on the value of a, and that no "correct" value 
of a could be specified so as to remove the bias [1^ . Since 
the time coordinate is an arbitrary label in the spacetime, 
we may impose the requirement that a viable measure 
prescription be invariant with respect to choosing even 
more general time gauges, such as 



Tim, 



(4) 



where T(0) > is an arbitrary function of the inflaton 
field (and possibly of other fields) , and the integration is 
performed along comovin g w orldlines x 



1,2,3 



const. 



The "spherical cutoff" [27[ and the "stationary mea- 
sure" [111 prescriptions were motivated by the need to 
remove the bias inherent in the equal-time cutoff. In par- 
ticular, the spherical cutoff selects as a compact subset V 
the interior of a large sphere drawn within the reheating 
surface R around a randomly chosen center. The spheri- 
cal cutoff is manifestly gauge-invariant since its construc- 
tion uses only the intrinsically defined 3-volume of the 
reheating surface rather than the spacetime coordinates 
(i, x). Some results were obtained in the spherical cutoff 
using numerical simulations (2^ . A disadvantage of the 
spherical cutoff is that its direct implementation requires 
one to perform costly numerical simulations of random- 
walk infiation on a spaceti me g rid, for instance, using the 
techniques of Refs. 0, [H, \EM- Instead, one would pre- 
fer to obtain a generally valid analytic formula for the 
probability distribution of cosmological observables. For 
instance, one could ask whether the results of the spher- 
ical cutoff depend in an essential way on the spherical 
shape of the region, on the position of the center of the 
sphere, and on the initial conditions. Satisfactory an- 
swers to these questions (in the negative) were obtained 



in Refs. [2^, l27[ in some tractable cases where results 
could be obtained analytically. However, it is difficult to 
analyze these questions in full generality since one lacks a 
general analytic formula for the probability distribution 
in the spherical cutoff. 

The RV measure is similar in spirit to the spherical 
cutoff because the RV cutoff uses only the intrinsic geo- 
metrical information defined by the reheating surface. It 
can be argued that the RV cutoff is "more natural" than 
other cutoffs in that it selects a finite portion V of the 
reheating surface without using artificial constant-time 
hypersurfaces, spheres, worldlines, or any other extrane- 
ous geometrical data. Instead, the selection of V in the 
RV cutoff is performed using a certain well-defined se- 
lection of subensemble in the probability space, which is 
determined by the stochastic evolution itself. 

The central concept in the RV cutoff is the "finitely 
produced volume." The basic idea is that there is always 
a nonzero probability that a given initial i/-region S does 
not give rise to an infinite reheating surface in its comov- 
ing future. For instance, it is possible that by a rare 
coincidence the infiaton field rolls towards reheating at 
approximately the same time everywhere in S. Moreover, 
there is a nonzero (if small) probability p{V)dV that the 
total volume Vol(i?) of the reheating surface R to the 
future of S belongs to a given interval [V, V -(- dV], 



p(V) = lim 



Prob{Vol(i?) e [V,V + dV]} 

dv ■ 



(5) 



I call p(V) the "finitely produced volume distribution." 
This distribution is nontrivial because the probability of 
the event Vol(i?) < oo is nonzero, if small, for any given 
(non-reheated) initial region S. The distribution p(V) is, 
by construction, normalized to that probability: 



p{V)dV = Prob{Vol(i?) < oo} < 1. 



(6) 



The RV cutoff consists of a selection of a certain en- 
semble of the histories that produce a total reheated 
volume equal to a given value V starting from an initial 
i/-rcgion. In the limit of large V, the ensemble con- 
sists of i7-regions that evolve "almost" to the regime of 
eternal infiation. Thus, heuristically one can expect that 
the ensemble -Ey provides a representative sample of the 
infinite reheating surface.'^ Given the ensemble E^, one 
can determine the volume-weighted probability distribu- 
tion p{Q\£v) of a cosmological parameter Q by ordinary 
sampling of the values of Q throughout the finite volume 
V. Finally, the probability distribution p((5) is defined as 
the limit of p(Q\£\;) at V — > oo, provided that the limit 
exists. 



^ Of course, this heuristic statement cannot be made rigorous since 
there exists no natural measure on the infinite reheating surface. 
We use this statement merely as an additional motivation for 
considering the RV measure. 
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To clarify the construction of the ensemble E^, it is 
helpful to begin by considering the distribution p{V) in 
a model that does not permit eternal inflation. In that 
case, the volume of the reheating surface is finite with 
probability 1, so the distribution p{V) is an ordinary 
probability distribution normalized to unity. In that con- 
text, the distribution p(V) was introduced in the recent 
work .14 .] where the authors considered a family of infla- 
tionary models parameterized by a number f2, such that 
eternal inflation is impossible in models where f2 > 1. It 
was then found by a direct calculation that all the mo- 
ments of the distribution p(V) diverge at the value = 1 
where the possibility of eternal inflation is first switched 
on. One can show that the finitely produced distribution 
p(V) for < 1 is again well-behaved and has finite mo- 
ments (see Sec HlDt . This FPRV distribution p(V) is the 
formal foundation of the RV cutoff. It is worth empha- 
sizing that the RV cutoff docs not regulate the volume 
of the reheating surface by modifying the dynamics of a 
given inflationary model and making eternal inflation im- 
possible. Rather, finite volumes V are generated by rare 
chance (i.e. within the ensemble Ey) through the un- 
modified dynamics of the model, directly in the regime 
of eternal inflation. 

Below I compute the distribution p{V) asymptotically 
for very large V in models of slow-roll inflation (Scc lIIDjl . 
Specifically, I will compute the distribution p(V;0o)j 
where (pQ is the (homogeneous) value of the inflaton field 
in the initial region S. To implement the RV cutoff ex- 
plicitly for predicting the distribution of a cosmological 
parameter Q, it is necessary to consider the joint finitely 
produced distribution p(V, Vg^; 0o, Qo) for the reheat- 
ing volume V{R) and the portion Vq^ of the reheating 
volume in which Q — Qn. (As before, (po and Qq are 
the values in the initial 7J- region.) If the distribution 
/9(V, Vq„; 00: Qo) is found, one can determine the mean 
volume (Vg^ly) while the total reheating volume V is 
held fixed. 



p(V;0o,Qo) 



(7) 



Then one computes the probability of finding the value of 
Q within the interval [Qr, Qr + dQ] at a random point 
in the volume V, 



pIQ^Qr-^v) 



V 



(8) 



The RV cutoff defines the probability distribution p{Q) 
for an observable Q as the limit of the distribution 
p{Q; V) at large V, 



p{Q) = lim 

V^oo 



V 



(9) 



One expects that this limit is independent of the ini- 
tial values (j)o,Qo because the large volume V is gener- 
ated by regions that spent a very long time in the self- 
reproduction regime and forgot the initial conditions. 



In Ref. [53| I derived equations from which the dis- 
tributions yo(V, Vq^; (/)o, Qo) and p{V;(j)o,Qo) can be in 
principle determined. However, a direct computation of 
the limit V — *■ oo (for instance, by a numerical method) 
will be cumbersome since the relevant probabilities are 
exponentially small in that limit. One of the main re- 
sults of the present article is an analytic evaluation of 
the limit V — > oo and a derivation of a more explicit for- 
mula, Eq. p7|) . for the distribution p{Q). The formula 
shows that the distribution p{Q) can be computed as a 
ground-state eigcnfunction of a certain modified Fokker- 
Planck equation. The explicit representation also proves 
that the limit ([9]) exists, is gauge- invariant, and is inde- 
pendent of the initial conditions (po and Qq. 

It was argued qualitatively in Ref. [53] that the RV 
measure does not suffer from the youngness paradox. In 
this article I demonstrate the absence of the youngness 
paradox in the RV measure by an explicit calculation. 
To this end, I will consider a toy model where every 
i?-region starts in the fluctuation-dominated (or "self- 
reproduction") regime with a constant expansion rate Hq 
and proceeds to reheating via two possible channels. The 
first channel consists of a short period dti of determinis- 
tic slow-roll inflation, yielding A'^i e-folds until reheating; 
the second channel has a different period St2 ^ Sti of de- 
terministic inflation, yielding N2 e-folds. (For simplicity, 
in this model one neglects fluctuations that may return 
the field from the slow-roll regime to the self-reproduction 
regime, and thus the time periods 6ti and St2 are sharply 
defined.) Thus there are two types of reheated regions 
corresponding to the two possible slow-roll channels. The 
task is to compute the relative volume- weighted probabil- 
ity P{2) /P{1) of regions of these types within the reheat- 
ing surface. (Essentially the same model was considered, 
e.g., in Refs. [H, HH, [ll [13 . See Fig. [H for a sketch of 
the potential l^((/))in this model.) 

This toy model serves as a litmus test of measure pre- 
scriptions. The "holographic" or "worldline" prescription 
yields P(2)/P(l) equal to the probability ratio of exiting 
through the two channels for a single comoving worldline. 
This probability ratio depends on the initial conditions. 
Thus, the worldline measure is (by design) blind to the 
volume growth during the slow-roll periods. On the other 
hand, the volume- weighted prescriptions of Refs. [2l|, [S^l 
both yield 



P(2) _ exp(3jV2) 
P(l) " exp(3Afi)' 



(10) 



rewarding the reheated iJ-regions that went through 
channel j by the additional volume factor exp(3Ar,). This 
ratio is now independent of the initial conditions. For 
comparison, an equal-time cutoff gives 



P(2) exp [3N2 - (3i/„ 



Ti - T2) St2] 



P(l) exp [3iVi - (3i7„,ax - Fi - Fa) Sti] ' 



(11) 



The overwhelming exponential dependence on dti and St2 
manifests the youngness paradox: Even a small difference 
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dt2—Sti in the duration of the slow-roll inflationary epoch 
leads to the exponential bias towards the "younger" uni- 
verses. The bias persists regardless of the choice of the 
time gauge [23!|, essentially because the presence of Sti 
and St2 in the ratio P(2)/P(l) cannot be eliminated by 
using a different time coordinate.* One expects that the 
RV measure will be free from this bias because the RV 
prescription does not involve the time coordinate t at all. 
Below (Sec. |IlF| I will show that the ratio P(2)/P(l) 
computed using the RV cutoff is indeed independent of 
the slow-roll durations Sti 2- The RV-regulated result 
[shown in Eq. ([55]) below] depends only on the gauge- 
invariant quantities such as A^i and A'2 and is, in general, 
different from Eq. PH)) . A calculation for an analogous 
landscape model was performed in Ref. [sl], yielding a 
result qualitatively similar to Eq. 

These calculations confirm that the RV measure has 
the desirable properties expected of a volume-based mea- 
sure: coordinate invariance, independence of initial con- 
ditions, and the absence of the youngness paradox. Thus 
the RV measure is a promising solution to the long- 
standing problem of obtaining probabilities in models of 
eternal inflation. Ultimately, the viability of the RV mea- 
sure proposal will depend on its performance in various 
example cases. In the calculations available so far, it is 
found that RV measure yields results that do not identi- 
cally coincide with the results of any other measure pro- 
posal. Hence, the RV measure is not equivalent to earlier 
proposals and needs to be studied in detail. 

As formulated here and in Ref. [s^, the RV measure 
prescription is directly applicable only to comparisons 
of reheating volumes, or in general of terminal states in 
the landscape (such as the anti-de Sitter bubbles). The 
RV proposal needs to be extended to predicting distribu- 
tions of properties not directly related to terminal states, 
such as the relative number of observations performed 
in different nonterminal bubbles. Then it will be possi- 
ble to investigate whether the RV measure suffers from 
the "Boltzmann brain" problem or from other difficulties 
encountered by some previous measure proposals. 

An extension of the RV measure to landscape scenar- 
ios can be achieved in several ways. For instance, one 
can consider the set of all possible future evolutions of a 
single nonterminal bubble and define the ensemble En of 
evolutions yielding a finite total number N of daughter 
bubbles (of all types). One can also consider the en- 
semble E'pf of evolutions yielding a finite total number 
N of observers in bubbles of all types. After computing 
the distribution of some desired quantity by counting the 
observations made within the finite set of N bubbles (or 
observers), the cutoff parameter N can be increased to 



* It should be noted that the youngness bias becomes very small, 
possibly even negligible, if one uses the number N of inflationary 
e-foldings as the time variable rather than the proper time t. I 
am grateful to A. Lindc and A. Vilenkin for bringing this to my 
attention. 



infinity. It remains to be seen whether the limit distribu- 
tions are different for differently defined ensembles, such 
as En and E'^, and if so, which definition is more suit- 
able. Future work will show whether some extension of 
the RV measure can provide a satisfactory answer to the 
problem of predictions in eternal inflation. 



II. OVERVIEW OF THE RESULTS 

In this section I describe the central results of this 
paper; in particular, I develop simplified mathematical 
procedures for practical calculations in the RV measure. 
For convenience of the reader, the results are stated here 
without proof, while the somewhat lengthy derivations 
are given in Sec. IIIII 



A. Preliminaries 

I consider a model of slow-roll inflation driven by an 
inflaton ch with the action 



-gd X. 



(12) 



In the semiclassical stochastic approach to inflation,^ the 
semiclassical dynamics of the field (f) averaged over an H- 
region is regarded as a superposition of a deterministic 
slow roll, 



(13) 



3H{(t>) AttG' 
and a random walk with root-mean-squared step size 

during time intervals St — H^^, where H((j)) is the func- 
tion defined by 



iJ(0) 



SttG 



Vi<t>). 



(15) 



A useful effective description of the evolution of the field 
at time scales St < H^^ can be given as 



(j){t + St) = (/)(<) + v{(l3)St + ^{t)y^2D{(t))St, (16) 
where ^(i) is a normalized "white noise" function, 

(0 = 0, mm) = sit-t'), (17) 



See Refs. [STlJSa. |59| for early works on the stochastic approach 
and Refs. [3,13 for pedagogical reviews. 
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which is approximately statistically independent between 
different 7?-regions. This stochastic process describes the 
evolution (j){t) and the accompanying cosmological expan- 
sion of space along a single comoving worldline. For sim- 
plicity, we assume that inflation ends in a given horizon- 
size region when = 0*, where 0* is a fixed value such 
that the relative change of H during one Hubble time 
5t = becomes of order 1, i.e. 









H 


4>='t>, 





^ 1. (18) 

4>=<t>, 



From the point of view of the stochastic approach, an 
inflationary model is fully specified by the kinetic coef- 
ficients I?(0), v{(f>), H{(f)). These coefficients are found 
from Eqs. p^ - p^ in models of canonical slow-roll in- 
fiation and by suitable analogues in other models. 

Dynamics of any fluctuating cosmological parameter Q 
is described in a similar way. One assumes that the value 
of Q is homogeneous in iJ-regions. The evolution of Q is 
described by an effective Langevin equation, 

Q{t + 6t) = Q{t) + vq{cp, Q)5t + ^Q{t)^2DQ{<i>, Q)St, 

(19) 

where the kinetic coefficients Dq and vq can be com- 
puted, similarly to D and v, from first principles. For 
simplicity we assume that the "noise variable" is in- 
dependent of the "noise" ^ used in Eq. (fTH]) . A corre- 
lated set of noise variables can be considered as well (see 
e.g. Ref. [H). 



where for brevity we dropped the subscript in (pg and 
Qo- This basic equation, first derived in Ref. [la], is 
of reaction-diffusion type and can be viewed as a non- 
linear modification of the Fokker-Planck equations used 
previously in the literature on the stochastic approach to 
infiation. 

While X{<I),Q) = 1 is always a solution of Eq. ([2T|l . 
it is not the correct one for the case of eternal infiation. 
A nontrivial solution, X{(j),Q) ^ 1, exists and has small 
values X{(j), Q) <gC 1 for (/>, Q away from the thermaliza- 
tion boundary. If the coefficients D/H and v/H happen 
to be Q-independent, the solution of Eq. (PTjl will be also 
independent of Q, i.e. X{(j),Q) = X{(j)), and thus deter- 
mined by a simpler equation obtained from Eq. pTjl by 
omitting derivatives with respect to Q, 

^^,00 + ^^,0 + 3X In X = 0. (23) 

It is easy to see that Eqs. (pij) and are manifestly 
gauge- invariant. Indeed, a change of time variable ac- 
cording to Eq. (HI) results in dividing the coefficients 
D, Dq,v,vq, H by the function T{(j}) [2J|, which leaves 
Eqs. (EB) and ([231) unmodified. 

Some approximate solutions of Eq. (|23p were given in 
Ref. [l^, where it was shown that X{(f>) is typically ex- 
ponentially small for in the inflationary regime. While 
small, X{(j)) is never zero; hence, there is a well-defined 
statistical ensemble of initial 7?-regions that have a finite 
total reheating volume in the future. The construction 
of the RV measure relies on this fact. 



B. Probability of finite inflation 

Let us consider an initial If-region S where the infiaton 
field as well as the parameter Q are homogeneous and 
have values (f> — (j)o and Q — Qo- For convenience we 
assume that reheating starts when = 0* and the Planck 
energy scales are reached at = 0pi independently of the 
value of Q. (If necessary, the field variables 0,(5 can be 
redefined to achieve this.) 

Although eternal inflation to the future of S is almost 
always the case, it is possible that reheating is reached at 
a finite time everywhere to the future of S, due to a rare 
fluctuation. In that event, the total reheating volume V 
to the future of S is finite. The (small) probability of 
that event, denoted by 

Prob(V < oo|0o,Qo) = ^(0o,Qo), (20) 

can be found as the solution of the following nonlinear 
equation, 

^X,^^ + ^X,QQ + ^ + ^XQ + 3X\nX = 0, 

(21) 

X(0pi,Q) = l, X(0*,Q) = 1, 

(22) 



C. Finitely produced volume 

In a scenario where eternal infiation is possible, we now 
consider the probability density p(V; 0o) of having a fi- 
nite total reheating volume V to the comoving future of 
an initial iJ-region with homogeneous value (j) — (fo- 
cusing attention at first on the case of infiation driven 
by a single scalar field). The distribution p(V; 0o) is nor- 
malized to the overall probability X(0o) of having a finite 
total reheating volume, 

/ p(V;0o)dV = X(0o). (24) 
Jo 

The distribution p(V; 0o) can be calculated by first deter- 
mining the generating function g{z; 0o), which is defined 
by 

poo 

g{z; 0o) EE (e-^^)^<^ ^ / e-^^p(V; 0o)dV. (25) 

This generating function is a solution of the nonlinear 
Fokker-Planck equation, 

Lg + 3glng = 0, (26) 
where the differential operator L is defined by 

L EE ^d^d^ + ^d^. (27) 
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In the case of several fields, say 0and Q, one needs to use 
the corresponding Fokker-Planck operator such as 



L 



QQ 



do dr. 



V4, 



dc 



The boundary conditions for Eq. (pS)) are 
g{z\(f),Q) = 1 for {0, Q} e Planck boundary, 



l{z;4>,Q) 



{0,Q} preheating boundary 



(28) 



(29) 
(30) 



Note that the parameter z enters the boundary condi- 
tions but is not explicitly involved in Eq. P^ . Also, the 
operator L and Eq. (|26p are manifestly gauge-invariant 
with respect to redefinitions of the form 

The generating function g plays a central role in the 
calculations of the RV cutoff. It will be shown below that 
the solution g{z; ip, Q) of Eq. needs to be obtained 
only at an appropriately determined negative value of z. 
This solution can be obtained by a numerical method or 
through an analytic approximation if available. 



D. Asymptotics of p(V; ^o) 

The finitely produced distribution p{V; (f) can be found 
through the inverse Laplace transform of the function 

9{z;(t)), 



1 /■'°° 
piV;(b) = — dz, 

27ri J-ioc 



(31) 



where the integration contour in the complex z plane 
can be chosen along the imaginary axis. The asymptotic 
behavior of p(V; (f)) at large V is determined by the right- 
most singularity of g{z] (j)) in the complex z plane. It 
turns out that the function g{z](j)) always has a singular- 
ity at a real, nonpositive 2; = of the type 

giz; (jj) = g{z,;4>) + a{^)^I^ + 0{z - z,), (32) 

where and (t((/)) are determined as follows. One con- 
siders the (z-dependent) linear operator 



L = i + 3(lng(z;</.) + l), 



(33) 



where L is the Fokker-Planck operator described above. 
For z > this operator is invertible in the space of func- 
tions satisfying zero boundary conditions. The value 
of z^ turns out to be the algebraically largest real number 
(in any case, z* < 0) such that there exists an eigenfunc- 

tion cr(0) of L with zero eigenvalue and zero boundary 
conditions, 



la{4>) = 0, 



(34) 



The specific normalization of the eigenfunction a{(\)) can 
be derived analytically but is unimportant for the present 
calculations. 



The singularity type shown in Eq. ([5^ determines the 
leading asymptotic of /o(V; 0) at V ^ cxd: 



P(V;0) 



2V^ 



z.V 



(35) 



The explicit form (|35p allows one to investigate the mo- 
ments of the distribution p(V;0). It is clear that all the 
moments are finite as long as z» < 0. However, if z* = 
all the moments diverge, namely for n > 1 we have 



(V"> 



p(V;0)V"dVcx / V"-3/2^V 



00. 



(36) 



The case z* = corresponds to the "transition point" 
analyzed in Ref. [14l |. corresponding to = 1 in their 



notation. This is the borderline case between the pres- 
ence and the absence of eternal inflation. The fact that 
z* = in the borderline case can be seen directly by not- 
ing that the Fokker-Planck operator L -I- 3 has in that 
case a zero eigenvalue, meaning that the 3-volume of 
equal-time surfaces does not expand with time (reheat- 
ing of some regions is perfectly compensated by inflation- 
ary expansion of other regions). In that case, the only 
solution g{z = 0; (/)) = A(0) of Eq. ([231) is A = 1 be- 
cause there are no eternally inflating comoving geodesies. 

Hence \Tig{z = 0, </>) = 0, and so the operator L is simply 

i = 1/ + 3. It follows that the operator L also has a zero 
eigenvalue at z = 0, and thus z z* = is the dom- 
inant singularity of g{z](f>). This argument reproduces 
and generalizes the results obtained in Ref. [ij where 
direct calculations of various moments of p{V; (f>) were 
performed for the case of the absence of eternal inflation. 

We note that the only necessary ingredients in the com- 
putation of (t((/)) is the knowledge of the singularity point 
z^ and the corresponding function g{z^] (/)), which is a so- 
lution of the nonlinear reaction-diffusion equation (|26p . 
Determining z* and g(z*; 0) in a given inflationary model 
does not require extensive numerical simulations. 



E. Distribution of a fluctuating field 

Above we denoted by Q a cosmological parameter that 
fluctuates during inflation but is in principle observable 
after reheating. One of the main questions to be an- 
swered using a multiverse measure is to derive the prob- 
ability distribution p{Q) for the values of Q observed in 
a "typical" place in the multiverse. I will now present a 
formula for the distribution p{Q) in the RV cutoff. This 
formula is signiflcantly more explicit and lends itself more 
easily to practical calculations than the expressions flrst 
shown in Ref. [ssj . 

As in the previous section, we assume that the dynam- 
ics of the inflaton field (f> and the parameter Q is de- 
scribed by a suitable Fokker-Planck operator L, e.g. of 
the form (pS)) . and that reheating occurs at = 0* 
independently of the value of Q. We then consider 
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Eq. for the function g{z;(l), Q) and the operator 

L = L + 3(lng + 1); we need to determine the value 

at which g{z;4>,Q) has a singularity. The operator L 
has an eigenfunction with zero eigenvalue for this value 
of z. This eigenfunction fa{z^;(f>,Q) needs to be deter- 
mined with zero boundary conditions (at reheating and 
Planck boundaries). Then the RV-regulated distribution 
of Q at reheating is 



p{Qii) — const 



dfo{z^;(f>, Q) D^^e 



-z,H- 



(37) 

where the normalization constant needs to be chosen such 
that / piQRjdQn = 1. The derivation of this result oc- 
cupies Sec. IIIIDI 

We note that /o is the eigenfunction /o of a gauge- 
invariant operator, and that the result in Eq. (1371) de- 
pends on the kinetic coefhcients only through the gauge- 
independent ratio D/H times the volume factor . 
The distribution p{Qr) is independent of the initial con- 
ditions, which is due to a specific asymptotic behavior of 
the finitely produced volume distributions, as shown in 
SecHlDl 



F. Toy model of inflation 

We now apply the RV cutoff to the toy model described 
at the end of Sec. HI We consider a model of inflation 
driven by a scalar field with a potential shown in Fig. [S] 
For the purposes of the present argument, we may as- 
sume that there is exactly zero "diffusion" in the deter- 
ministic regimes ipi^'' < cj) < and 02 < < 0*^\ while 
the range </>! < < 02 is sufficiently wide to allow for 
eternal self-reproduction. Thus there are two slow-roll 
channels that produce respectively A^i and N2 e-folds 
of slow-roll inflation after exiting the self-reproduction 
regime. Since the self-reproduction range generates arbi- 
trarily large volumes of space that enter both the slow-roll 
channels, the total reheating volume going through each 
channel is inflnite. We apply the RV cutoff to the prob- 
lem of computing the regularized ratio of the reheating 
volumes in regions of types 1 and 2. 

In this toy model it is possible to obtain the results 
of the RV cutoff using analytic approximations. The 
required calculations are somewhat lengthy and can be 
found in Sec. IIIIEI The result for a generic case where 
one of the slow-roll channels has many more e-folds than 
the other (say, N2 ^ Ni) can be written as 



m 
p(i) 



0(1) 



H 



-3/a(2) 



) exp [3iV2] 



ij-3(0(i))cxp[3iVi] 



exp [3iVi2] , (38) 



where we have defined 
^2 
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V2Hi 



Hi ^ ^-^Vo. (39) 




Figure 2: A model potential with a flat self-reproduction 
regime 01 < < 02 and deterministic slow-roll regimes 
0!^'' < < 01 and 02 < < 0!^^ producing A''i and A'^2 
inflationary e-folds respectively. In the interval 0i < < 02 
the potential V{(j}) is assumed to be constant, V{(j>) ~ Vo- 



The pre-exponential factor 0(1) can be computed nu- 
merically, as outlined in Sec. IIII El 

We note that the ratio ([38|l is gauge-invariant and does 
not involve any spacetime coordinates. This result can be 
interpreted as the ratio of volumes e'^^^ and e^^^ gained 
during the slow-roll regime in the two channels multiplied 
by a correction factor e'^^^^ . The dimensionless number 
A^i2 can be suggestively interpreted (up to the factor -\/2) 
as the mean number of "steps" of size S(f) ^ required 
for a random walk to reach the boundary of the flat region 
[01 J 02] starting from the middle point 0o = ^ (0i +02)- 
Since each of the "steps" of the random walk takes a Hub- 
ble time Hq^ and corresponds to one e- folding of infla- 
tion, the volume factor gained during such a traversal 
will be e'^^^^. Note that the correction factor increases 
the probability of channel 2 that was already the domi- 
nant one due to the larger volume factor e'^^^ ^ e'^^^. 
Depending on the model, this factor may be a significant 
modification of the ratio PH)) obtained in previously used 
volume-based measures. 



III. DERIVATIONS 
A. Positive solutions of nonlinear equations 

It is not easy to demonstrate directly the existence of 
nontrivial solutions of reaction-diffusion equations such 
as Eq. (pij) . However, there is a connection between so- 
lutions of such nonlinear equations and solutions of the 
linearized equations. Rigorous results are available in the 
mathematical literature on nonlinear functional analysis 
and bifurcation theory. 

Heuristically, consider a solution of Eq. that is 
approximately ^^^(0) ~ 1. The equation can be linearized 
in the neighborhood ofX«lasX==I — x{4>) ^nd yields 
the Fokker-Planck (FP) equation 



L + 3 



X = 0, L = ^d^^ + ^d4,. (40) 



The FP operator L + 3 is adjoint to the operator 
+ 3 



P^d^^[^p]-d^{^P) + 3P, (41) 
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which enters the FP equation for the 3-volume distri- 
bution P{4>, t) in the e-folding time parameterization. If 
eternal inflation is allowed in a given model, the operator 
+ 3 has a positive eigenvalue. The largest eigenvalue 
of that operator is zero in the borderline case when eter- 
nal inflation is just about to set in. The spectrum of the 
operator L + 3 is the same as that of the adjoint oper- 
ator U + 3. Hence, in the borderline case the largest 
eigenvalue of the operator L + 3 will be zero, and there 
will exist a nontrivial, everywhere nonnegative solution 
X of Eq. (pOj) . Thus, heuristically one can expect that 
a nontrivial solution X((j)) ^ 1 will exist away from the 
borderline case, i.e. when the operator L+3 has a positive 
eigenvalue. 

Following the approach of Ref . , one can imagine a 
family of inflationary models parameterized by a label Q, 
such that eternal inflation is allowed when < 1 . Then 
Eq. (f23|) will have only the trivial solution, X{4i) = 1, for 
57 > 1. The case f2 = 1 where eternal inflation is on the 
borderline of existence is the bifurcation point for the so- 
lutions of Eq. At the bifurcation point, a nontrivial 
solution X(</)) ^ 1 appears, branching off from the trivial 
solution. A rigorous theory of bifurcation can be devel- 
oped using methods of nonlinear functional analysis (see 
e.g. chapter 9 of the book (S^)- In particular, it can be 
shown that a nontrivial solution of a nonlinear equation, 
such as Eq. (j23p , exists if and only if the dominant eigen- 
value of the linearized operator L + 3 with zero boundary 
conditions is positive. 

There remains a technical difference between the eigen- 
value problem for the operator L + 3 with zero boundary 
conditions and with the "no-diffusion" boundary condi- 
tions normally used in the stochastic approach. 



d 



[D(0)P(0)] = 0. 



(42) 



It was demonstrated in Ref. [l5j that the eigenvalue of 
L + 3 with the boundary conditions (|42)) is positive if a 
nontrivial solution of Eq. (j23p exists. In principle, the 
eigenvalue of L -I- 3 with zero boundary conditions is not 
the same as the eigenvalue of the same operator with the 
boundary conditions ((42|) . One can have a borderline 
case when one of these two eigenvalues is positive while 
the other is negative. In this case, the two criteria for the 
presence of eternal inflation (based on the positivity of 
the two different eigenvalues) will disagree. However, the 
alternative boundary conditions are imposed at reheat- 
ing, i.e. in the regime of very small fluctuations where the 
value of the eigenfunction P(0) is exponentially small 
compared with its values in the fluctuation-dominated 
range of </>. Hence, the difference between the two eigen- 
values is always exponentially small (it is suppressed at 
least by the factor e~^^, where is the number of e-folds 
in the deterministic slow-roll regime before reheating). 
Therefore, we may interpret the discrepancy as a limita- 
tion inherent in the stochastic approach to inflation. In 
other words, one cannot use the stochastic approach to 



establish the presence of eternal inflation more precisely 
than with the accuracy e""^^. Barring an extremely fine- 
tuned borderline case, this accuracy is perfectly adequate 
for establishing the presence or absence of eternal infla- 
tion. 

The main nonlinear equation in the calculations of 
the RV cutoff is Eq. for the generating function 

g{z;(f)). That equation differs from Eq. (j23p mainly by 
the presence of the parameter z in the boundary con- 
ditions. Therefore, solutions of Eq. (|26p may exist for 
some values of z but not for other values. Note that 
g{z = 0;(/)) = A(0); hence, nontrivial solutions g{z;(j)) 
exist for z — Q under the same conditions as nontriv- 
ial solutions X{(l)) ^ 1 of Eq. P5|) . While it is certain 
that solutions g{z\ (p) exist for z > 0, there may be val- 
ues z < for which no real- valued solutions g(z;(j)) exist 
at all. However, the calculations in the RV cutoff re- 
quire only to compute g{z^,; cj)) for a certain value z.^, < 0, 
which is the algebraically largest value z where g{z; (jj) 
has a singularity in the z plane. The structure of that 
singularity will be investigated in detail below, and it will 
be shown that g{z^,; (f) is finite while dg/dz oc (z— z*)^"'^/^ 
diverges at z = z*. Hence, the solution g{z;(f)) remains 
well-defined at least for all real z in the interval [z*, -|-oo]. 
It follows that g{z; (f>) may be obtained e.g. by a numeri- 
cal solution of a well-conditioned problem with z = z^+e, 
where e > is a small real constant. 



B. Nonlinear Fokker-Planck equations 

In this section I derive Eq. (I^S]) . closely following the 
derivation of Eq. (pij) in Ref. pL5j] . 

We begin by considering the case when inflation is 
driven by a single scalar field (p, such that reheating 
is reached at (j) = Let /o(V;0o) be the probability 
density of obtaining the finite reheated volume V. We 
will derive an equation for a generating function of the 
distribution of volume, rather than an equation directly 
for piV'jipo)- Since the volume V is by definition non- 
negative, it is convenient to define a generating function 
g{z; (f>o) through the expectation value of the expression 
exp(— zV), where z > is the formal parameter of the 
generating function. 



5(z;0o)^(e-^^) 



V<oo 



e-'^p{V;Po)dV. (43) 



Note that for any z such that Rez > the integral in 
Eq. converges, and the events with V = -|-oo are 
automatically excluded from consideration. However, we 
use the subscript "V < oo" to indicate explicitly that 
the statistical average is performed over a subset of all 
events. The distribution p(V; 4>o) is not normalized to 
unity; instead, the normalization is given by Eq. 

The parameter z has the dimension of inverse 3- 
volume. Physically, this is the 3-volume measured along 
the reheating surface and hence is defined in a gauge- 
invariant manner. If desired for technical reasons, the 
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variable z can be made dimensionless by a constant 
rescaling. 

The generating function g{z\ (f>) has the following mul- 
tiplicative property: For two statistically independent re- 
gions that have initial values (p — (j)i and (j) — (j)2 respec- 
tively, the sum of the (finitely produced) reheating vol- 
umes Vi -f V2 is distributed with the generating function 

(44) 

We now consider an _ff-region at some time having 
an arbitrary value 0(<) not yet in the reheating regime. 
Suppose that the finitely produced volume distribution 
for this i7-rcgion has the generating function g[z] </>). Af- 
ter time 5t the initial iJ-region grows to iV = e^^^*" sta- 
tistically independent, "daughter" i?-regions. The value 
of 4) in the fc-th daughter region (fc = 1, A^) is found 
from Eq. 



(45) 



where the "noise" variables (fc = 1, N) are statisti- 
cally independent because they describe the fluctuations 
of 4> in causally disconnected 7J-regions. The finitely pro- 
duced volume distribution for the fc-th daughter region 
has the generating function g{z\4)k)- The combined re- 
heating volume of the A'" daughter regions must be dis- 
tributed with the same generating function as reheating 
volume of the original iJ-region. Hence, by the multi- 
plicative property we obtain 



N 



g 



(46) 



k=l 



We can average both sides of this equation over the noise 
variables to get 



N 



(47) 



\k=l 



Since all the ^fc are independent, the average splits into 
a product of A^ identical factors, 



7(2:; 0) 



g(z;0 + ^;(0)<5^-t-^/2^e) 



TV 



(48) 



The derivation now proceeds as in Ref. [if 
compute, to first order in 5t, 



We first 



g{z- + v{(t))5t + v/2i?(0)O) = 5 + (^^5,0 + Dg^^^) St. 

(49) 

Substituting A^ — e^^^^ and taking the logarithmic 
derivative of both sides of Eq. (|48|) with respect to 6t 
at St = 0, we then obtain 



0= Aln5(z;0) 



3H\ng 



^'g,0 + Dg,< 
g 



(50) 



The equation for g(z; 0) follows, 

Dg,4,4, + w.g,0 + 3Hg \ng = 0. (51) 
This is formally the same as Eq. (PTjl . However, the 



boundary conditions for Eq. (j51[) are different. The con- 
dition at the end-of-inflation boundary = 0,^ is 



g(z;0*) = e 



.) 



(52) 



because an i?-region starting with = 0* immediately 
reheats and produces the reheating volume -ff~'^(0*). 
The condition at Planck boundary 0pi (if present), or 
other boundary where the effective field theory breaks 
down, is "absorbing," i.e. regions that reach = 0pi do 
not generate any reheating volume: 



5(2; 0: 



PI) 



1. 



(53) 



The variable z enters Eq. (|5ip as a parameter and only 
through the boundary conditions. At z = the solution 
isg(O;0)=X(0). 

A fully analogous derivation can be given for the gen- 
erating function g{z] 0o, Qo) in the case when additional 
fluctuating fields, denoted by Q, are present. The gener- 
ating function g(z; 0o, Qo) is defined by 



7(z;0o,(9o) = / e 
/o 



-zV 



p(V;0o,Qo)dV, (54) 



where p(V; 0o, Qo) is the probability density for achieving 
a total reheating volume V in the future of an ff-region 
with initial values 0o,Qo of the fields. In the general 
case, the fiuctuations of the fields 0, Q can be described 
by the Langevin equations 



(t>{t + St) = (t>{t) + v^St + ^0 y/2D,p^St + y/2DpQSt, 

(55) 



Q{t + 6t) = Q{t) + vqSt + S.4,yj2D^QSt + ^Q^2DQQSt, 

(56) 

where the "diffusion" coefficients D^^f,, D^q, and Dqq 
have been introduced, as well as the "slow roll" veloci- 
ties Vff, and VQ and the "noise" variables and ^q. The 
resulting equation for (7(z;0o,(3o) is (dropping the sub- 
script 0) 

Lg + Sgliig = 0, (57) 
where the differential operator L is defined by 



H 



H 



H 



H 



(58) 

The ratios D^^/H, etc., are manifestly gauge-invariant 
with respect to time parameter changes of the form Q . 

Performing a redefinition of the fields if needed, one 
may assume that reheating is reached when = 0* in- 
dependently of the value of Q. Then the boundary con- 
ditions for Eq. ([57|l at the reheating boundary can be 
written as 



g{z\ 



,Q) 



zH' 



(59) 



The Planck boundary still has the boundary condition 
5(z;0Pi) = 1. 
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C. Singularities of g{z) 

For simplicity we now focus attention on the case of 
single- field inflation; the generating function g{z;(j)) then 
depends on the initial value of the inflaton field (j). The 
corresponding analysis for multiple fields is carried out 
as a straightforward generalization. 

By definition, g{z; (j)) is an integral of a probability dis- 
tribution p(V; (j)) times e~^^. It follows that g{z; cf)) is an- 
alytic in z and has no singularities for Re z > 0. Then the 
probability distribution p(V; (f) can be recovered from the 
generating function g(z;(j)) through the inverse Laplace 
transform, 



1 r°° 

piV;(b)^— dz, 
27ri J_ioo 



(60) 



where the integration contour in the complex z plane can 
be chosen along the imaginary axis because all the sin- 
gularities of g{z] (f>) are to the left of that axis. The RV 
cutoff procedure depends on the limit of p(V; (/)) and re- 
lated distributions at V ^ oo. The asymptotic behavior 
at V — *■ cxD is determined by the type and the location 
of the right-most singularity of g{z; (j)) in the half-plane 
Re z < 0. For instance, if z = z* is such a singularity, the 
asymptotic is p(V; 4>) cx exp[— z*V]. The prefactor in this 
expression needs to be determined; for this, a detailed 
analysis of the singularities of g(z; (j)) will be carried out. 

It is important to verify that the singularities of g(z; 0) 
are (/)-independent. We first show that solutions of 
Eq. ((57)) cannot diverge at finite values of ip. If that were 
the case and say g{z]4)) — > oo as ^ (/>i, the function 
Inlng as well as derivatives g^^ and g^^^ would diverge 
as well. Then 



lim dfjj [In In g\ = lim 



d,j,g 
g^'ag 



oo. 



(61) 



It follows that the term glng is negligible near (p — (pi 
in Eq. ([57]) compared with the term d^g and hence also 
with the term dtpd^g. In a very small neighborhood of 

= </>! , the operator L can be approximated by a linear 
operator Li with constant coefficients, such as 



L w Li = Aid4,d4, + Bid^ 



(62) 



Since at least one of the coefficients Ai, Bi is nonzero at 
(j) — (j)i/it follows that g{z\ (j)) is approximately a solution 
of the linear equation Lig — near (j) = (pi. However, 
solutions of linear equations cannot diverge at finite val- 
ues of the argument. Hence, the function g(z; (p) cannot 
diverge at a finite value of (p. 

The only remaining possibility is that the function 
g{z;<p) has singular points z — z* such that g{z^,;(p) re- 
mains finite while dg/dz, or a higher-order derivative, 
diverges at z = z*. We will now investigate such diver- 
gences and show that g{z;(p) has a leading singularity of 
the form 



where z* is a (/>-independent location of the singularity 
such that z* < 0, while the function (7(0) is yet to be 
determined. 

Denoting temporarily gi(z;0) = dg/dz^ we find a lin- 
ear equation for gi , 



L5i-f3(lng+l)5i ==0, 



(64) 



with inhomogeneous boundary conditions 

gi((/),) = -i/-3(0,)e-^^"'(**), 5i(M = 0. (65) 

The solution .gi(z; (p) of this linear problem can be found 
using a standard method involving the Green's function. 
The problem with inhomogeneous boundary conditions is 
equivalent to the problem with zero boundary conditions 
but with an inhomogeneous equation. To be definite, let 
us consider the operator L of the form used in Eq. ([5T|) . 



H{cP) 



dsdd, 



(66) 



Then Eqs. ([64|) -([65 |) are equivalent to the inhomogeneous 
problem with zero boundary conditions, 

Lgi + 3 (Ing + 1) 51 = DH-^e-'"'" 5' [cP - 0,), (67) 
5i(z;0*) = 5ri(z;(/)pi) = 0. (68) 

The solution of this inhomogeneous equation exists as 
long as the linear operator L -I- 3(ln g + 1) does not have 
a zero eigenfunction with zero boundary conditions. 

Note that the operator L -I- 3 (In g -I- 1) is explicitly z- 
dependent through the coefhcient g{z\ (p). Note also that 
g{z;(p) ^ by definition (|^ for values of z such that 
the integral in Eq. converges; hence \ng is finite for 
those z. Let us denote by G(z; (p, (p') the Green's function 
of that operator with zero boundary conditions. 



LG + 3 (In .g(z, 0) -f 1) G = Si<p - cp') 
G{z-4>.A')^G{z-^^,, 



0. 



(69) 
(70) 



This Green's function is well-defined for values of z such 
that L + g{z\ (p) + 1) is invertible. For these z we may 
express the solution gi{z;(j)) of Eqs. explicitly 
through the Green's function as 



5i(z;0) = - 



zH- 



dG{z-cP,<P') 



(71) 



9{z;(p') = 9{z*;(p') + (Ti4>)y/z - z^ + 0{z - z*), (63) 



Hence, for these z the function gi{z\ (p) = dg/dz remains 
finite at every value of <p. A similar argument shows 
that all higher-order derivatives d^g/dz"^ remain finite 
at every (j) for these z. Therefore, the singularities of 
g{z; (p) can occur only at certain (/)-independent points 
z — z^ , z — z^ , etc. 

Since the generating function g(z; (p) is nonsingular for 
all complex z with Rez > 0, it is assured that gi{z;(p) 
and G(z; (p, cp') exist for such z. However, there will be 
values of z for which the operator L -t- 3(ln g -I- 1) has a 
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zero eigenfunction with zero boundary conditions, so the 
Green's function G is undefined. Denote by such a 
value with the algebraically largest real part; we already 
know that Re z.^ < in any case. Let us now show that 
the function gi{z;(j)) actually diverges when z z^. In 
other words, lim^^^^ gii^', (j)) — c>o for every value of (j). 

To show this, we need to use the decomposition of the 
Green's function in the eigenfunctions of the operator 
L + 3(ln.g + l), 

G(z; <t>, 0') = £ ^/„(0)/*(0'), (72) 

where fn{z', (j)) are the (appropriately normalized) eigen- 
functions with eigenvalues A„(z) and zero boundary con- 
ditions. 



L + 3(lng(z; cj>) + 1)J /„(z; 0) = A„(z)/„(z; 4>), (73) 

fn{z\ 0) = for = = (/)pi. 

(74) 



The decomposition (f72|) is possible as long as the oper- 
ator L is self-adjoint with an appropriate choice of the 
scalar product in the space of functions /(0). The scalar 
product can be chosen in the following way, 



(/l,/2 



/l(0)/2*(<^)M(0)d0, 



(75) 



where M{(j)) is a weighting function. One can attempt to 
determine M{(j)) such that the operator L is self- adjoint. 



(/l,i/2) = (L/1,/2) 



(76) 



In single-field models of inflation where the operator L 
has the form (|66t , it is always possible to choose M{(p) 
appropriately [241 • However, in multi-field models this 
is not necessarily possible.^ One can show that in stan- 
dard slow- roll models with K fields (j)K and kinetic 
coefficients 



D, 



ij3 



r2 ' 



1 dH 



there exists a suitable choice of M{(j)), namely 



M{(bi,...,4 
such that the operator 



, ttG 
K) = exp 



ttG 

1p 



(77) 
(78) 

(79) 



I am grateful to D. Podolsky for pointing this out to me. The 
hermiticity of operators of diffusion type in the context of eternal 
inflation was briefly discussed in Ref. |6ll . 



is self-adjoint in the space of functions /(</>) with zero 
boundary conditions and the scalar product (j75p . How- 
ever, the operator L may be non-self-adjoint in more gen- 
eral inflationary models where the kinetic coefficients are 
given by different expressions. We omit the formulation 
of precise conditions for self-adjointness of L because this 
property is not central to the present investigation. In 
non-self-adjoint cases a decomposition similar to Eq. ([7^ 
needs to be performed using the left and the right eigen- 
functions of the non-self- adjoint operator L + 3(ln g +1). 
One expects that such a decomposition will still be pos- 
sible because (heuristically) the nondiagonalizable oper- 
ators are a set of measure zero among all operators. The 
requisite left and right eigenfunctions can be obtained 
numerically. We leave the detailed investigation of those 
cases for future work. Presently, let us focus on the case 
when the decomposition of the form fTS]) holds, with 
appropriately chosen scalar product and the normalized 
eigenfunctions 



{fni: fri) — 



(80) 



The eigenfunctions fmiz; (j)) can be obtained e.g. numer- 
ically by solving the boundary value problem ([75)l - d7i)) . 

In the limit z ^ z*, one of the eigenvalues A„ ap- 
proaches zero. Since linear operators such as L al- 
ways have a spectrum bounded from above (2^ , we may 
renumber the eigenvalues A„ such that Aq is the largest 
one. Then we define z* as the value with the (alge- 
braically) largest real part, such that Ao(z*) = 0. By 
construction, for all z with Re z > Re z* all the eigenval- 
ues A„ are negative. Note that the (algebraically) largest 
eigenvalue Ao(z) is always nondegenerate, and the corre- 
sponding eigenfunction /o(z;0) can be chosen real and 
positive for all 0, except at the boundaries (/> = 0* and 
(j) — (jj-pY where / satisfies the zero boundary conditions. 

For z near z*, only the nondegenerate eigenvalue Aq 
will be near zero, so the decomposition ([7^ of the Green's 
function will be dominated by the term I/Aq. Hence, we 
can use Eqs. ([7T|) and ([7^ to determine the function 
gi{z;(f)) approximately as 



51(2;'/') 



faiz;4>) 9/o(z;0^) 
Ao(z) d(j) 



D 

IF' 



00 as z 



n) 



because 



It follows that indeed gi{z;(/) 
Ao(z) ^ 0. 

This detailed investigation allows us now to determine 
the behavior of g{z; 4>) at the leading singularity z — z^. 
We will consider the function g{z\4>) for z near z* and 
show that the singularity indeed has the structure 

We have already shown that the function g itself does 
not diverge at z = z^ but its derivative gi = dg/dz does. 
Hence, the function g{z^](j)) is continuous, and the dif- 
ference g{z^,;(f)) — g{z]<j)) is small for z « z*, so that we 
have the expansion 

5g{z-4,)=g{z-cp)-g{z,-4>) (82) 
^gi{z-(^){z-z,) + 0[{z,-zf]. (83) 
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(Note that we are using the finite value gi{z](l)) rather 
than the divergent value gi (z* ; 0) in the above equa- 
tion.) On the other hand, we have the explicit repre- 
sentation ([5T|) . Let us examine the values of Xo{z) for 
z « z*. At z — we have Ao(z*) = 0, so the (small) 
value Ao(z) for z « can be found using standard per- 
turbation theory for linear operators. If we denote the 
change in the operator L by 



dL = 3(In ,g(z; (/)) - lng(z*; (/<)) 
we can write, to first order, 

Ao(z)« (/o,<5L/o\ = / |/o(<^)| 



_ 3'5g(z;0) 
3(2:*; 0) '' 

2 3(5.9 (z; 0) 



(84) 



9{z*]4>) 



Now, Eqs. (EI]) and ^ yield 



Sg{z;(p) _ fo{z;<P) 
z — z^, Ao(z) 



dch m 



d(j3. (85) 



(86) 



Integrating the above equation in <j) with the prefactor 

3 



i/o(0)r 



(87) 



and using Eq. ([85]) . we obtain a closed equation for Ao(z) 
in which terms of order (z — z.^)^ have been omitted. 



Ao(^) 



Ao(z) 



86 H^' 



I r ,2 SSg{z;(j3) 
\Jo\ —7 ttJo 



It follows that Ao(z) cx y/ z — z^, and .gi(z;0) cx 
(z — z*) , confirming the leading asymptotic of the 
form 

Let us also obtain a more explicit form of the singular- 
ity structure of ^(z;^). We can rewrite Eq. ([88]) as 



Ao(z) sa aoVz - z* + 0(z - z*). 



(89) 



where co is a constant that may be obtained explicitly. 
Then Eq. jH]) yields 



gi{z;<p) « 

Y z z* 



(90) 



with a different constant ai. Finally, we can integrate 
this in z and obtain 



g{z;^) = g{z,;^) + 2aifo{z;^)^I~i:+0{z-z,). (91) 
We may rewrite this by substituting z z* into /o(z; 0), 



g{z; (j)) = g{z.,; 0) + a{(f>)^/ z - z, + 0{z - z,), (92) 
a{cj)) = 2aifoiz,;c^). (93) 



The result is now explicitly of the form (|63p . It will turn 
out that the normalization constant 2ai cancels in the fi- 
nal results. So in a practical calculation the eigenfunction 
/o(z*;0) may be determined with an arbitrary normal- 
ization. 

As a side note, let us remark that the argument given 
above will apply also to other singular points z^ ^ z^ as 
long as the eigenvalue Afe(z) of the operator i7-|-3(ln g-f 1) 
is nondegenerate when it vanishes at z = z^. If the 
relevant eigenvalue becomes degenerate, the singularity 
structure will not be of the form z — z'^ but rather 
(z — z^)* with some other power < s < 1. 

Now we are ready to obtain the asymptotic form of the 
distribution p(V; </)) for V — > 00. We deform the integra- 
tion contour in the inverse Laplace transform (j60p such 
that it passes near the real axis around z = z* . Then we 
use Eq. ([92)) for g{z; (f>) and obtain the leading asymptotic 



p(V;0) 



1 
27ri 
1 



a(0) 
a(0)V 



\/ z — z^ e dz 



(94) 



The subdominant terms come from the higher-order 
terms in the expansion in Eq. (j92p and are of the order 
V^^ times the leading term shown in Eq. (|94p . 

Finally, we show that z* must be real-valued and that 
there are no other singularities z^ with Re z^ — Re z* . 
This is so because the integral 



51 (z*;, 



p(V;0)e- 



VdV 



00 



(95) 



will definitely diverge for purely real z* if it diverges for 
a nonreal value z^ = z.^ -f iA. If, on the other hand, the 
integral (I95|) diverges for a real z* , it will converge for any 
nonreal z^ = z* -l-iA with real A because the function 
p(V; 4>) has the large-V asymptotic of the form ([M)) and 
the oscillations of exp(i74V) will make the integral 
convergent. 



D. FPRV distribution of a field Q 

In this section we follow the notation of Ref. [H^ . 

Consider a fluctuating fleld Q such that the Fokker- 
Planck operator L is of the form ([55]) . We are inter- 
ested in the portion Vg^ of the total reheated volume 
V where the field Q has a value within a given interval 
[Qr,Qr + dQ]. We denote by p(V,Vqh;(/)o,(9o) the joint 
probability distribution of the volumes V and Vq^ for ini- 
tial if-regions with initial values cf) — (pQ and Q = Qq. 
The generating function g(z, q; cf), Q) corresponding to 
that distribution is defined by 



g{z,q;(f>,Q) 



-'^v-<iVQp{y,VQ-^,Q)dVdVQ. (96) 



Since this generating function satisfies the same mul- 
tiplicative property (|44p as the generating function 
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g[z; (j), Q), we may repeat the derivation of Eq. ([57)1 with- 
out modifications for the function g{z, q; (j), Q). Hence, 
g{z, q; (j), Q) is the sohition of the same equation as 
g{z;(f),Q). The only difference is the boundary condi- 
tions at reheating, which are given not by Eq. ((59l) but 
by 

g{z,q;(j)„Q)^exp[-{z + q6QQ„)H-^{(j)„Q)] , (97) 

where (with a shght abuse of notation) Sqqj^ is the indi- 
cator function of the interval [Qr, Qr + dQ], i.e. 



Sqq^ = e{Q - QR)e{QR + dQ- Q). 



(98) 



We employ this "finite" version of the J-function only be- 
cause we cannot use a standard Dirac (5-function under 
the exponential. This slight technical inconvenience will 
disappear shortly. 

The solution for the hmction g{z, q; (j), Q) may be ob- 
tained in principle and will provide complete information 
about the distribution of possible values of the volume Vq 
together with the total reheating volume V to the future 
of an initial i/-region. In the context of the RV cutoff, 
one is interested in the event when V is finite and very 
large. Then one expects that Vq also becomes typically 
very large while the ratio Vq/V remains roughly con- 
stant. In other words, one expects that the distribution 
of Vq is sharply peaked around a mean value (Vq), and 
that the limit (Vq) /V is well-defined at V ^ oo. The 
value of that limit is the only information we need for 
calculations in the RV cutoff. Therefore, we do not need 
to compute the entire distribution p(V, Vg; (j), Q) but only 
the mean value ( VqI^) at fixed V. 

Let us therefore define the generating function of the 
mean value ( VQ|y) as follows. 



/j(z;<^,g)^(VQ,e-^^) 



V<^ Qq^^^l 



O;0,Q). 



(99) 

(The dependence on the fixed value of Qr is kept im- 
plicit in the function /i(z; 0, Q) in order to make the no- 
tation less cumbersome.) The differential equation and 
the boundary conditions for h{z\ (j), Q) follow straightfor- 
wardly by taking the derivative 9^ at g = of Eqs. ((57|) 
and ([97|l . It is clear from the definition of g that g{z, q ~ 
0; 0, Q) = g{z; 4>, Q). Hence we obtain 



L/i + 3(ln5(z;0,Q) + l)h^O, 
h{z;(f)^,Q) 



h{z;4)pi,Q) = 0. 



(100) 

-SQQn, (101) 
(102) 



Note that it is the generating function g, not g, that 
appears as a coefficient in Eq. (jlOOp . 

Since the "finite" (5-function <5qqh now enters only lin- 
early rather than under an exponential, we may replace 
^QQr by the ordinary Dirac (5-function 5{Q — Qr). To 
maintain consistency, we need to divide h by dQ, which 
corresponds to computing the probability density of the 



reheated volume with Q = Qr. This probability density 
is precisely the goal of the present calculation. 

The RV-regularized probability density for values of Q 
is defined as the limit 



p{QFt) = „lim 



lim 



(Vq«Iv) 



V/V<oo 



Vp{V-^,Q) 



plot 



h{z\ (j), Q)dz 



v^oo V Je'^^g{z;(j),Q)dz 



(103) 



To compute this limit, we need to consider the asymp- 
totic behavior of {^QrIv)^^^ at V ^ oo. This behav- 
ior is determined by the leading singularity of the func- 
tion h{z; (j), Q) in the complex z plane. The arguments of 
Sec. nil Cl applv also to h{z] (/), Q) and show that h cannot 
have a </>- or Q-dependent singularity in the z plane. 

Moreover, h{z] (j), Q) has precisely the same singular 
points, in particular z = z*, as the basic generating 
function g{z; </>, Q) of the reheating volume. Indeed, the 
function h{z; (j>, Q) can be expressed through the Green's 
function G{z; (j), Q, (j>' , Q') of the operator L + 3(lng + 1), 
similarly to the function 51(2:; 0) considered in Sec. IIII Cl 
For z z^,, this operator is invertible on the space of func- 
tions /(0, Q) satisfying zero boundary conditions. Hence, 
h{z; (j), Q) is nonsingular at z 7^ and becomes singular 
precisely at z = z*. 

Let us now obtain an explicit form of h{z; (p, Q) near 
the singular point z ~ z^. We assume again the eigen- 
function decomposition of the Green's function (with the 
same caveats as in Sec. IIII C"|) . 

00 

G(z; 0, Q, 0', Q') - VTT/"('^; Q'), 

(104) 

where /„ are appropriately normalized eigenfunctions of 
the z-dependent operator L -\- 'iihi g -\- 1) with eigenvalues 
A„(z). The eigenfunctions /„ must satisfy zero boundary 
conditions at reheating and Planck boundaries. Similarly 
to the way we derived Eq. ((8T|). we obtain the explicit 
solution 



/i(z;0,Q) = ^ 



/„(z;0, Q) 
A,i(z) 



dfn D^e' 



zH- 



(105) 

The value of h{z; (j), Q) for z sa z.^ is dominated by the 
contribution of the large factor l/Ao(z) cx (z — z*) , 
so the leading term is 



h{z; cf), Q) 



fo{z^;(t),Q) 
Ao(z) 



dfo D^se-'*"' 



ij4 



(106) 

The asymptotic behavior of the mean value ( "I^Qr lv)v<oo 
is determined by the singularity of h{z; (f), Q) at z = z*. 
As before, we may deform the integration contour to pass 
near the real axis around z = z*. We can then express 
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the large-V asymptotic of ( "^'qr lv)v<oo follows, 



k{z*](j), Q) 



,z,V 



9/0 DA,d,e~ 



.,Qr 
(107) 



where ctq is the constant defined by Eq. 

We now complete the analytic evaluation of the 
limit (|103p . Since the denominator of Eq. (|103p has the 
large-V asymptotics of the form 



V / g{z- 0, Q)e^'^dV cx /o(z,; 0, Q)V—-e 



(108) 



where /o is the same eigenfunction, the dependence on 
and Q identically cancels in the limit (|103p . Hence, that 
limit is independent of the initial values cf) and Q but 
is a function only of Q/j, on which h{z] 0, Q) implicitly 
depends. Using this fact,, we can significantly simplify 
the rest of the calculation. It is not necessary to compute 
the denominator of Eq. (jlOSp explicitly. The distribution 
of the values of Q at </> = 0^ is simply proportional to 
the Q /f-dependent part of Eq. (|107p ; the denominator of 
Eq. (jl03[) serves merely to normalize that distribution. 
Hence, the RV cutoff yields 



const 



d(h iJ4 



(109) 

where the normalization constant needs to be chosen such 
that / p{QR)dQR = 1- This is the final analytic formula 
for the RV cutoff applied to the distribution of Q at re- 
heating. The the value z.^,, and the corresponding solution 
g{z^;(j), Q) of Eq. ([26)) . and the eigenfunction fo{z^;(t>, Q) 
need to be obtained numerically unless an analytic solu- 
tion is possible. 

Let us comment on the presence of the factor D^^ in 
the formula (|109p . The "diffusion" coefficient is eval- 
uated at the reheating boundary and is thus small since 
the fluctuation amplitude at reheating is (in slow-roll in- 
flationary models) 
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-5 



(110) 



Nevertheless it is not possible to set Z?^^ = directly in 
Eq. (|109p . This is so because the existence of the Green's 
function of the Fokker-Planck operator such as L depends 
on the fact that L is a second-order differential operator 
of elhptic type. If one sets D^^ = near the reheating 
boundary, the operator L becomes first-order in </) at that 
boundary. Then one needs to use a different formula than 
Eq. (j67p for reducing an equation with inhomogeneous 
boundary conditions to an inhomogeneous equation with 
zero boundary conditions. Accordingly, one cannot use 



formulas such as Eq. ([7T|) for the solutions. Alternative 
ways of solving the relevant equations in that case will 
be used in Sec. IIIIEI 



E. Calculations for an inflationary model 

In this section we perform explicit calculations of RV 
cutoff for a model of slow-roll inflation driven by a scalar 
field with a potential shown in Fig. [2l The kinetic coeffi- 
cients D{(p) and v{(p) are such that D{(j)) — Dq, v{(p) — 0, 
and H{(p) = Ho in the flat region (pi < (f> < (f>2, where 
the constants Dq and Hq are 



-Vn. 



Dn = 



-"0 
87r2' 



(111) 



In the slow-roll regions 0i < < (f>i^^ and <^2 < (f> < (0\ 
the coefficient 1^(0) is set equal to zero, while «(</>) 7^ 
and i?(0) is not constant any more. The number of e- 
folds in the two slow-roll "shoulders" can be computed by 
the standard formula, 



-d<p = -4^G / —c 
V Ja,, H' 



J = 1,2. 



(112) 

The first step of the calculation is to determine the 
singular point z = of solutions g{z;(j>) of Eq. ([5T|) . 
We expect 2* to be real and negative. The boundary 
conditions for g(z](j)) are 



exp 



zH- 



In each of the two deterministic regions, 
and (f)2 < (p < 4>* 1 Eq- l|5T|) becomes 



V 

H 



3g\ng = 0, 



with the general solution 



g{z; 4>) = exp 




Cexp -3 / — 



(113) 



< < 01 



(114) 



(115) 



where C is an integration constant. Since the equation 
is first-order within the deterministic regions, the solu- 
tions are fixed by the boundary condition (|113p in the 
respective region. 



g{z; 0) = exp 



-zH 




)exp -3/ — 



(116) 



We may therefore compute the values of 5(2; 0) at the 
boundaries 0i_2 of the self-reproduction region as 

g[z- 01,2) = exp Uzi?-3(^(i,2)) (3^^^^)1 ^ (^^7) 
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Now we need to solve Eq. (IFT|) with these boundary con- 
ditions in the region c/fi < < 02 ■ The equation has then 
the form 



(118) 



Exact solutions of Eq. (|118p were studied in Ref. [l5| , to 
which the reader is referred for more details. It is easy 
to show that Eq. (jllSp is formally equivalent to a one- 
dimensional motion of a particle with coordinate g{(t>) in 
a potential U{g), 



U{g) = %g' (21n.9-l), 



(119) 



while (j) plays the role of time. A solution g{z; (p) with 
boundary conditions (jll7[) corresponds to a trajectory 
that starts at the given value g{z]4>i) with the initial 
velocity chosen such that the motion takes precisely the 
specified time interval 4)2 — 0i and reaches g(z\(\)2). For 
z < the boundary conditions specify g{z;(l)i,2) > 1, 
i.e. the trajectory begins and ends to the right of the 
minimum of the potential (sec Fig. ^ . Since the system 
is conservative, there is a constant of motion E (the "en- 
ergy") such that E = U{go) at the highest point of the 
trajectory go where the "kinetic energy" vanishes. The 
solution g{z; (p) can be written implicitly as one of the 
two alternative formulas, 



± 



dg 



^2E{z) - 2U{g) 



(120) 



valid in appropriate intervals 4>i < 4> < po and 4>o < (f) < 
4>2 respectively, where (po is the value of cj) corresponding 
to the turning point g^ — g{z; po). The value E = E{z) 
in Eq. (|120p must be chosen such that the total "time" is 







•'go 


J ga 



dg 



y/2E ~ 2U{g) 



(121) 



This condition together with E = U (go) implicitly deter- 
mine the values E = E{z) and go = go{z)- 

The singularity z = of the solution g{z] (p) is found 
by using the condition dg/dz oo. Differentiating 
Eq. (|120p with respect to z and substituting Eq. (|117p 
for g{z] (pi.2), we obtain the condition 



dg 



1 



dz ^2E{z) - 2U{g) ^2E{z) - 2U{g{z-p,a)) 



Ig{z;4>) [2E{z) - 2U{g)r^ 
It follows that dg/dz ^ 00 when 

E{z)^U{giz;Pi^2)) 



= 0. 



(122) 



(123) 



This condition is interpreted in the language of Fig. [3] as 
follows. As the value of z becomes more negative, the 



Uig) 
90 



0.5 



9{z;4>i) g{z;p2) 
1/ 1^5^ 



E 



6tt' 
' Hi, 




Figure 3: Potential U((j) given by Eq. (|119|l can be used to 
interpret solutions g(z\ (p) as mechanical motion in "time" p at 
constant total energy E. The potential vanishes at <; = and 
g = e^^^ and has a minimum at g = 1. The value go is the 
turning point where U{go) = E. A trajectory corresponding 
to z < will begin and end with g > 1, i.e. to the right 
of the minimum of the potential. Solutions cease to exist 
when z < z^; the solution g{z,;p), shown by the thin line 
with arrows, corresponds to a trajectory that starts at rest 
(as demonstrated in the text). The energy of this trajectory 
is i5 ~ 0, and so the value g{z,;p2) is close to e^^^ while 
g{z,;pi) is close to 1. 



initial and the final values of g given by Eq. (|117p both 
grow. The last available trajectory starts from rest at 
4> = (p2 and at the value of g such that U{g) — E. 

To obtain a specific result, let us assume that 4>2 — 4>i is 
sufficiently large to provide self-reproduction [p2 — 0i 3> 
Ho) and that the number of e-folds in channel 1 is smaller 
than that in channel 2, 



H- 



6l^Vxp(3Afi) < i?-^(?!)?Vxp(37V2). (124) 



Then the value g{z; P12) will grow faster than g{z; pii) as 
z becomes more negative. It follows that g{z] 4>2) will 
reach the singular point first. Since the "time" 02 — (t>i 
is large, the constant E will be close to so that the 
trajectory spends a long "time" near g = 0. Then the 
value g{zs,;(p2) will be close to e^^^. Hence the value of 
z^ is approximately 

w -ii/3(^(2)) gj.p (_37V2) . (125) 

For this value of z.^,, the starting point of the trajectory 
will be 



g(z*; pi) ss exp 



i exp (37Vi - 3A^2) 



(126) 
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Hence, the solution 17(2:*; (f>) at the singular point z = 
can be visualized as the thin line in Fig. [31 starting ap- 
proximately at g{z^,;(f>i) = 1 and finishing at 5(0*; 02 ) ~ 

An approximate expression for ; </>) can be obtained 
by setting £' « in Eq. ()120p : then the integral can be 
evaluated analytically. In the range (/>o < < (/'2 we 
obtain 



dg 



so the solution is 



Hp 
Vl2^ 



v/l - 21ng 







02 

(127) 



g{z^\(t)) ~ exp 



&0 < < 



In the range 0i < 
approximation 



(128) 

< 00 we obtain within the same 



.g(z*; 0) « exp 



1 6tt^ ( 



Ho 



These approximations are valid for within the indicated 
ranges and away from the turning point 0o. The value 
of 00 can be estimated by requiring that the value of 
(7(z»;0o) obtained from Eq. (|128p be equal to that ob- 
tained from Eq. (fT29| . This yields 



Ha 



(130) 



We note that the value g{z^; 0o) can be obtained some- 
what more precisely by approximating the solution 
g{za.;4>) in a narrow interval near = 0o by a function 
of the form exp [A + B{(j) — 0o)^] and matching both the 
values and the derivatives of g{z^,;(f)) to the approxima- 
tions (|128p and (|129p at some intermediate points strad- 
dling (f) — (j)Q. In this way, a uniform analytic approxima- 
tion for 5(2;*; 0) can be obtained. However, the accuracy 
of the approximations (|128p and (|129p is sufficient for the 
present purposes. 

Having obtained adequate analytic approximations for 
z* and (7(2:*, 0), we can now proceed to the calculation of 
the mean volumes (Vi,2|y) of regions reheated through 
channels 1 and 2 respectively, conditioned on the event 
that the total volume of all reheated regions is V. We use 
the formalism developed in Sec. lIIIDl where the variable 
Q now takes only the discrete values 1 and 2, so instead 
let us denote that value by j. The relevant generating 
fimction hj{z](j)) is defined by 

/i,(z;0)^(V,e-^^)^<^, J = 1,2, (131) 



and is found as the solution of Eq. (|100p . which now takes 
the form 



v{4>) 



+ 3(ln.9(z;0) + l) 



hjiz;(t)) = 0, 
(132) 



with boundary conditions imposed at the reheating 
boundaries, 

/ii(z;0W) = H-3(0«)e-^-^(^l^'), /.i(z;0l^)) = 0; 

(133) 

/i2(z;0l^)) = 0, /^2(z;0l^)) = i7-^(0?))e-^^-^(^^^'). 

(134) 

In the present toy model the diffusion coefficient is set to 
zero at reheating, so the formalism developed in Sec. lIII Dl 
needs to be modified. We will first solve Eq. (|132p ana- 
lytically in the no-diffusion intervals of and obtain the 
boundary conditions for h at the boundaries of the self- 
reproduction regime [0i,02] where D{(f>) 7^ 0. Then the 
methods of Sec. IIII Dl will be applicable to the boundary 
value problem for the interval [0i, 02]- 

Implementing this idea in the first no-diffusion region 
^ ^ ^ fj^Q solution pi6p for g{z;(j)) and 



< < 



reduce Eq. to 



(129) d,h,+ 



V 



1 - zH' 



exp 



jj 
— c 

(1) V 



h, = 0. 



(135) 



This equation is easily integrated together with the 
boundary conditions (I133p - (|134p and yields the values 



of hj at 



/ij(z;0i) =5,i77-^(0l'^)exp 



3A^i - zH 



(136) 

Similarly we can determine the values hj (z; 02) . Since the 
value z = z* is important for the present calculation, we 
now find the values of hj sd z = z^ using the assumption 
N2 > iVi and the estimate (|125l) . 



hj{z^;(j)i) w 6i 



exp [3Ni 



i/3( 



i,j = l,2. (137) 



We have thus reduced the problem of determining 
hj(z;(f>) to the boundary- value problem for the interval 
[01 1 4'2] where the methods of Sec. IIII Dl apply but the 
boundary conditions are given by Eq. p37p . 

The next step, according to Sec. IIIIDl is to compute 
the cigcnfmiction /o(z*; 0) of the operator 

1 = :^5^9^ + 3(ln5(z,;0) + l) (138) 
^0 



such that 



i/o-0; /o(z,;0i,2) = 0. (139) 



As we have shown, this eigenfunction with eigenvalue 
exists precisely at z = z*. Once this eigenfunction is 
computed, the ratio of the RV-regulated mean volumes in 
channels 1 and 2 will be expressed through the derivatives 
of /o at the endpoints and through the modified boundary 
conditions p37p as follows. 



Pm ^ h2{z,;<f>2) |g0/o(z,;02)| 
P(l) ft,i(z*;0i) 90/o(z*;0i) ' 



(140) 
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Figure 4: Sketch of the eigenfunction fo{z,;(l>) (dashed hne) 
interpreted as the wavefunction of a stationary state with zero 
energy in the potential V{(l)) (solid line). Due to exponential 
suppression by the potential barrier, the amplitude of /o in 
the right region is exponentially larger than that in the left 
region. 



a potential barrier separating two classically allowed re- 
gions near (j) = and (p = 4>2- Since the "potential well" 
a,t (j) — (j)2 is deeper, the ground state is approximately 
the ground state of that one well, with an exponentially 
small amplitude of being near cj) = <f>i. The shape of the 
eigenfunction is sketched in Fig.[4l The exponential sup- 
pression of the amplitude near (j) — (j)i can be found using 
the WKB approximation, which yields 



l^0/o(^*;02)| 



= A21 cxp 



(143) 



where (f)i^2 are the turning points such that V{(f>i^2) — 0. 
The pre-exponential factor A21 is of order 1 and can, 
in principle, be obtained from a more detailed matching 
of the WKB-approximated solution across the barrier to 
the solutions in the "classically allowed" regions, or by 
determining the solution /o(^*;0) numerically. However, 
we will omit this calculation since the main result will 
consist of an exponentially large factor. That factor can 
be estimated using Eqs. (fTM]) . (fT29l) . and (fUTjl as 



The absolute value is taken to compensate for the neg- 
ative sign of the derivative d^fo at the right bound- 
ary point (assuming that fo > everywhere). Since 
^1,2 (-z*; '/'1, 2) are already known, it remains to derive an 
estimate for fo{z^,;<p). 

The eigenvalue equation L/q = formally resembles 
a one-dimensional Schrodinger equation with the coordi- 
nate (j) and the "potential" 



127r2 
^0 



(lng(z,;0)-M). 



(141) 



The eigenfunction /o(2;»;(/)) is then interpreted as the 
"wavefunction" of a stationary state with zero energy 
and zero boundary conditions at = 01.2- According to 
Eqs. (|128p and (|129p . the function V{(f>) has a maximum 
at « 00 (see Fig. \^ , while its values at the endpoints 
are 



V{<l>i 



T^(02 



IStt^ 



(142) 



Using the terminology of quantum mechanics, there is 



V{(l))d(t> « 2 



00 



127r2 , 



3n' 



H, 







Hence, the ratio (|140p is simplified to 



(144) 



P(2) . i/-3(0?h e3^^+5 



P(l) 



^21- 



H- 



■ exp 



(02 - (pl) 



(145) 



This is the main result quoted above in Eq. ([38 
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